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Abstract 

We revisit four generations within the context of supersymmetry. We compute the per- 
turbativity limits for the fourth generation Yukawa couplings and show that if the masses 
of the fourth generation lie within reasonable limits of their present experimental lower 
bounds, it is possible to have perturbativity only up to scales around 1000 TeV. Such low 
scales are ideally suited to incorporate gauge mediated supersymmetry breaking, where the 
mediation scale can be as low as 10-20 TeV. The minimal messenger model, however, is 
highly constrained. While lack of electroweak symmetry breaking rules out a large part of 
the parameter space, a small region exists, where the fourth generation stau is tachyonic. 
General gauge mediation with its broader set of boundary conditions is better suited to 
accommodate the fourth generation. 



1 Introduction 



Interest in a sequential fourth generation of chiral fermions has waxed and waned over the 
past decades with the changing status of constraints implied by various precision measure- 
ments in the flavour and the gauge sector. The recent observation of single-top events at 
the Tevatron [1,2] allowed a direct and clean determination of the CKM matrix element 
\V t b\ = 0.91 ± 0.11 (stat+syst) ± 0.07 (theory) [2]. This is in good agreement with the 
Standard Model prediction of |V#,| ~ 1, but falls short of excluding extensions of the the- 
ory by another chiral generation of quarks and leptons, as has been emphasised in recent 
publications [3 — 8] . As a matter of fact, the mixing between the third and a hypothetical 
fourth familjo can be as large as the mixing between the first two generations in the Stan- 
dard Model and yet be compatible with all experimental data including direct searches for 
new quarks and leptons, electroweak precision measurements, and flavour changing neutral 
currents [4-6,8,9]. 

Currently, we do not have a theoretical understanding of the number of generations, and 
a priori there is no reason why there should not be another one. A fourth generation 
would have profound implications for particle physics phenomenolog}{§. For one thing, 
it can ease the tension between the LEP bound on the Higgs mass and the electroweak 
precision measurements [4,6,8] as we will elaborate on in Section [2j Secondly, the Yukawa 
couplings of the fourth generation fermions, heavier than those of the first three, will have 
important implications for the perturbativity of the theory at the high scale. This can 
have nontrivial effects, for example, on the Higgs mass bounds obtained in the SM by 
demanding that the Landau pole lie above the Planck scale [11,12]. As a result of this 
possible effect of the fourth generation on the consistency of the theory up to high scale, 
unified theories with four generations have also received special attention, both with su- 
persymmetry (SUSY) [13-20] and without it [21]. The fourth generation can also account 
for the extra CP violation needed for electroweak baryogenesis to work [22-24]. Other 
ideas that have been explored in connection with four generations include: Extra dimen- 
sions [25-27], technicolour [28-30] and electroweak symmetry breaking [31-33], radiative 
mass generation [34,35], bounds from cosmology [36], neutrino physics [37-39], and finally 
string theory [40]. 

Recent analyses of the fourth generation have mainly focused on the the non-supersymmetric 
case. In this, we extend the minimal supersymmetric standard model (MSSM) by one chiral 
generation and explore the implications for supersymmetry and supersymmetry breaking. 
For clarity, we will call the MSSM with three and four generations MSSM3 and MSSM4, 
respectively. In Section [21 we review the constraints on the masses of the fourth gener- 
ation coming from experiment, precision electroweak data, and flavour changing neutral 
currents. In Section [31 we introduce our notation and map out the parameter space of the 
MSSM4 where the theory remains perturbative up to some assumed unification scale of 

■""We will denote the fourth generation quarks and leptons by b' , t', v' t . For enhanced readability, in 
graphs we may use the alternate notation i4, 64, t4, v t ±. 
2 For a review, see Rcf. [10]. 
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the order 2 x 10 16 GeV. We will find that this puts severe restrictions on m t >, my, m T >, 
tan p, and show that the current experimental bounds and perturbative unification are 
mutually exclusive. Depending on the masses and tan/3, the theory becomes strongly- 
coupled around 10-1000 TeV. To illustrate the qualitative differences, we will present a toy 
mSUGRA modeH where we have chosen the quark and lepton masses to be equal to their 
third generation counterparts. In Section HI motivated by the low "perturbativity" scale, 
we explore this issue in the context of gauge mediated supersymmetry breaking (GMSB). 
Minimal GMSB, however, suffers from tachyons in the spectrum, and hence we are led 
to generalise our GMSB set-up in the quest for realistic models. Finally, in Section [5] we 
summarise our results and outline directions for future work. 



2 Limits from Experiment, Electroweak Precision Data, 
and FCNC 

The current experimental limits quoted by the PDG [41] at 95% CL are : 

m t , > 256 GeV, my > 128 GeV, m T , > 100.8 GeV, m K > 45 GeV. (1) 

The bounds on m t > and my assume that the predominant decay mode is to a W boson 
and another quark [42,43], which we expect to be true for a sequential fourth generationE). 
Along the same lines, for the limit on m T i one assumes that r' decays to a W and a stable 
z/ [46], and the limits for a stable heavy neutral particle give the lower bound on m v > [47]. 

Electroweak precision measurements further constrain the allowed mass ranges. For no 
mixing between the third and fourth generation, x 2 is minimised for \m t > — my \ ~ 45 — 75 
GeV [4], and for m t > = 300 GeV with my subject to this constraint, the mixing can be 
as large as sin 6*34 = 0.35 [6]. The precision data excludes larger mixing between the third 
and fourth families [6] that is otherwise allowed from FCNC constraints [5]. 

Fig. 10.4 of Ref. [41] by the PDG shows the constraints from electroweak precision mea- 
surements. The values S = —0.04 ± 0.09 and T = 0.02 ± 0.09 obtained from the fit are 
in best agreement with a Higgs mass of m# = 117 GeV. For higher Higgs masses, the 
90% CL contours move towards smaller S and larger T values. From this it is clear that a 
fourth generation may ease the tension between the LEP bound on the Higgs mass and the 
electroweak precision data, since its contribution to T is typically positive [4,6]. A fourth 
generation can give a negative contribution to 5*, but the mass splitting of the quarks is 
constrained by the T parameter and leads to a positive correction to 5* that has to be kept 
small in order to stay in the 90% CL ellipse. 



Note that the limits quoted in Eq. (1) are at 95% CL. In the absence of direct availability 



of these bounds at a higher level of confidence, we approximate them by subtracting 20% 



3 Hcre and in the following, wc will use the word mSUGRA, where we should more correctly call it the 
constrained minimal supersymmetric standard model (CMSSM). 

4 One can find a higher lower bound on mi,/ by adding model-specific assumptions [44,45]. 
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off the respective limits at 95% CL. Furthermore, the exclusion limits denote the pole 
masses, whereas in our calculations, we need to use the running masses. We will account 
for this difference by taking yet another 5% off the masses for QCD corrections. Note that 
this is a conservative estimate as the SUSY threshold corrections can induce an additional 
difference up to 20% between the pole and the running masses [48]. 

In our analysis, we will thus be working with two sets of mass limits, namely those at 95% 
CL as well as the weaker lower bounds obtained with the above-mentioned prescription. 
Both sets are subject to the constraint from electroweak precision measurements that the 
mass splitting in the same SU(2) multiplet be not greater than ~ 75 GeV. Thus we consider 
in our analysis the following sets of fourth generation fermion masses: 

m t > = 256 GeV, m v = 181 GeV, m T > = 100.8 GeV (2) 

m t > = 192 GeV, m v = 117 GeV, m T > = 75 GeV (3) 
We will comment on the mass of the fourth generation neutrino in Section [3j 

3 Perturbativity and Four Generations 

One of the main constraints in considering models with four generations is the perturba- 
tivity of the Yukawa couplings. With the masses of the fourth generation expected to be 
typically larger than the already known third generation masses, the scale up to which the 
theory remains perturbative is a major concern. It is expected that supersymmetry would 
soften the running of the Yukawa couplings and enable the theory to be valid up to much 
higher energy scales. It is well known that the Yukawa couplings of the heavy fermions 
flow under the renormalisation group to a infrared (quasi-)fixed point. The proximity of 
the observed top mass to this limit (see e.g. Ref. [49]), in fact also indicates that a fourth 
generation of fermions, heavier than the top quark, can have important implications for 
the high scale perturbativity of the theory. 

In the present section, we discuss this issue in detail and present the results of our com- 
putations. As we will show, these high scales will extend only up to ~ 1000 TeV. Let us 
also note that in supersymmetric theories the scale up to which perturbativity is preserved 
has implications for (a) gauge coupling unification and (b) the scale of supersymmetry 
breaking, though a priori they are two independent sectors of the theory. 

The generalisation of the MSSM to four generations is straight-forward. We will denote the 
fourth generation as the primed generation: (f, b', t', u' t ) and denote this model as MSSM4. 
The standard MSSM picture with three generations will from now on be mentioned as 
MSSM3. The MSSM4 superpotential takes the form : 

W = u c Y U QH 2 + d c YdQtfi + e c Y e LH 1 + nH x H 2 , (4) 
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where Y are 4x4 matrices in generation space. The Yukawa couplings of the fourth 
generation, defined in terms of their masses, are given as 

m t 'V2 myy/2 m T iy/2 

h t < = — — j , h v = , h T > = -, (5) 

vsmp v cos p v cos p 

where v = 246 GeV stands for the Higgs vacuum expectation value. 

In addition, the fourth generation neutrino v T < should also attain a mass m v i > 45 GeV. 
There are several ways of generating the fourth generation neutrino mass term , which may 
either be of the Dirac type or a mixture of Dirac plus Majorana type, leading to a Majorana 
mass for the v T < and this can make the analysis highly model dependent. In the present 
case, we have not considered the effects of a Dirac Yukawa coupling for the neutrino in the 
renormalisation group equations (RGE). It should be noted that a large neutrino Dirac 
Yukawa coupling can, in principle, affect strongly the evolution of the Yukawa coupling of 
the t' but, the effect would be minimal as long as the said neutrino Dirac Yukawa coupling 
is small. A more detailed analysis of the various possible neutrino Yukawa couplings and 
their impact on perturbativity will be presented elsewhere. 

Firstly, let us note that the larger masses of the down type fourth generation fermions 
mean that requiring that the Yukawa couplings be perturbative at the weak scale puts 
upper bounds on tan j3, stronger than the ones present in the MSSM3. The strongest 
limitH comes from hy. Imposing that h v ~ Air, we have : 

i 

tan/? < (2n (v/myf - l) 2 (6) 



For the lower limit of my quoted in Eq. (2) on the previous page this already sets a limit 



tan/? < 4.7. The bound from m T i is much weaker, tan/3 ~ 38. Even for the values of 
tan/3 obeying this limit, the Yukawa couplings are not expected to be perturbative all the 
way up to the GUT scale, where the gauge coupling unification happens. 

It should be noted that in the presence of four generations, gauge coupling unification at 
the 1-loop level still takes place at the scale Mgut ~ 2 x 10 16 GeV, though the value 
of a GUT itself changes. We refer to Appendix I A. 1 1 for the relevant renormalisation group 
equations. At the 2-loop level, the renormalisation group running of gauge couplings will 
involve Yukawa couplings and thus if the Yukawa couplings become non-perturbative at 
scales much lower than M GUT they could render the same to the gauge couplings. In 
the non-supersymmetric Standard Model with four generations, it has been known that 
perturbative unification of the gauge couplings is possible [21]. However, the values for 
the fourth generation masses used in the analysis of Ref. [21] are ruled by the recent 
experimental results (see also Ref. [12]). 

In the following, we will derive upper limits on the fourth generation fermion masses by 
requiring that the Yukawa couplings of the fourth generation be perturbative all the way 
up to the GUT scaled 



5 Lower limits on mt> on the other hand put lower limits on tan/3 of limited consequence. 
6 Related analysis can be found in Refs. [19,20]. 
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For this purpose, we use the two loop RGE equations listed in Appendix IA. 21 for uq number 
of generations. Needless to say, hq = 3 for MSSM3 and uq = 4 for MSSM4. We solve 
these equations numerically for a given set of masses at the weak scale and check whether 
the corresponding Yukawa couplings remain perturbative at the high scale. 

We have developed a variant of the popular supersymmetric spectrum calculator, SOFTSUSY 
[50], called INDISOFT. It contains significant modifications including the ability to handle 
three or four generations. Some technical aspects are presented in the Appendix [Bj 

We first vary m t i and my while keeping m T i fixed and later vary my and m T > while keeping 
m t i fixed. From the RGE in Appendix IA. 21 we see that the evolution of t' is independent of 
the t' mass at 1-loop order. With this in mind, for the first analysis we take the r' mass to 
be negligible. For numerical purposes we set it equal to the r mass, i.e. m T > = m T . Before 
presenting the numerical results one final comment is in order. The Yukawa couplings of 
the heavy fields (e.g. t', t) flow towards an infrared fixed point at the weak scale. Using 
this one can derive an upper bound on these masses of the heavy fermions analytically. 
The analytical results for the t' are presented in Appendix IA.3I 



In the left column of Fig. 1 on the following page, we plot the regions in the my — m t 



plane where the Yukawa couplings remain perturbative all the way up to the GUT scale 
using 2-loop MSSM4 RGE. We have chosen tan/5 = 3. In the first row, the condition 
that only the top-prime Yukawa coupling remains perturbative all the way is plotted. The 
perturbativity limit is taken to be ~ 3.54. All the regions are colour-coded with the 
values that the Yukawa couplings take at the high scale. The legend for the colour-coding 
is shown on the right side of each plot. If the Yukawa coupling attains a value beyond the 
upper limit of 3.54, the point is flagged as non-perturbative and is denoted in black. Thus, 
in each of these plots regions in black are ruled out by the perturbativity limit. From this 
first plot of the left column, we see that a large region in m t i opens up as the b' mass 
increases. For my > 70 GeV, m t / > 150 also can be made valid. However, although the 
t' Yukawa coupling is perturbative in these regions, the other two Yukawa couplings are 
not, as is evident from the plots in the next two rows. In the second row, left column, we 
exhibit the same plots, however requiring that only hy remains perturbative, and in the 
third row, left column, requiring that only h T > remains perturbative. As we see from the 
plots, although the t' Yukawa coupling is perturbative, the h' and r' Yukawa couplings are 
no longer perturbative in these regions. The last row shows plots where all the constraints 
are put together. Here we see that for a negligible r' mass (~ 1.75 GeV), t' and b' masses 
are constrained to be: 

0<m b , <70GeV, < m t / < 160 GeV, m T i = 1.78 GeV, tan/5 = 3. (7) 

It is instructive to compare the above results with those of the MSSM3. In MSSM3, as we 
know, the Yukawa couplings do remain perturbative all the way up to the GUT scale for 
the known masses of the top and bottom quarks and the tau lepton, the only exception 
being in regions where tan/3 is very small or very large. In the right panel of Fig. 1, we 



present the analogous analysis for the case of the MSSM3. We have fixed the r mass at 
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Figure 1: The first three rows of the left and right panels show ht>, hy, h T > in MSSM4, and h t , hb, 
h T in MSSM3, respectively, at 2-loop and Mx = 2.3 x 10 16 GeV, plotted in the mt'-mb' and mt-mb 
plane for tan/3 = 3 and m T i = m T = 1.78 GeV. For the MSSM4 case this experimentally excluded value 
is chosen for illustrative purposes only. The black regions indicate where the Yukawa couplings become 
non-perturbative. The last row shows the regions wffere all Yukawa couplings are perturbative. 



its experimental value but varied m t and m&. We have further fixed tan/3 to be 3 as in 
the case of four generations. The results in the case of MSSM3 are strikingly similar to 
those obtained for the four generational MSSM4. In fact we can read off from the plot the 
regions in m t — m& plane where all the three Yukawa couplings remain perturbative: 

0<m b <60GeV, 0<m t <180GeV, tan/3 = 3. (8) 

The reason for these similar results is the way Landau poles appear in the Yukawa couplings. 
The evolution of the Yukawa couplings depends only on themselves and the gauge couplings. 
As it is clear from the RGEs (see Appendix IA.3j) . the evolution of the gauge couplings is 
very similar to that of the three generation case, except for the /3-functions b{ which change 
the slope of the h(t). The change in the hi for the three vs. four generations is not very 
large to induce large changes in the upper limits. The difference is within a few tens of 
GeV. 




(a) Constraints in the TOb'-mf plane from the per- (b) Constraints in the mb>-m T i plane from the 
turbativity of h t ' for fixed values of m T > = 100.8 perturbativity of the max (ht> , hy , h T >) , for fixed 
GeV and tan/3 = 3. values of mt> = 150 GeV and tan/3 = 3. 

Figure 2: MSSM with 4 generations at 2-loop. The black region is where the respective Yukawa couplings 
become non-perturbative below the unification scale Mx = 2.3 x 10 16 GeV. 



The results of Eq. (7) on page 6 are valid only when m T > = m T . Increasing the value of m T > 



to the lower limit of Eq. (2) on page 4 significantly modifies these results. We find that as 
expected mp which is less dependent on the t' mass remains perturbative for a similarly 
large region of the parameter space as shown in Fig. ^|(a) However, for the values, 



0<m*<0 GeV, < m t , < 160 GeV, 



100.8 GeV, tan/3 = 3, (9) 



hy is not perturbative for any my. In fact, h T > itself is no longer perturbative. 

We now proceed to keeping mp fixed while varying my and m T i. We have chosen two 
values of mp, the one given by the approximate perturbative upper limit ~ 150 GeV 
Eq. (7) on page 6) and the other the experimental lower limit of 256 GeV. In Fig. 2||(b) 
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we have plotted the regions in the my —m T i plane in which the theory remains perturbative 
all the way up to the GUT scale. From the figure we can read off the valid mass ranges as: 



0<m 6 , <60GeV, 0<m r , <50GeV, tan/? = 3, m t , = 150 GeV (10) 



0<m b ,<5QGeV, < m T > < 53 GeV, tan (3 = 3, 



m v 



256 GeV 



WW 



Clearly these mass ranges are excluded by the constraints of Eq. (1) on page 3 or even by 



the weaker constraints in Eq. (2) and Eq. (3) on pageHJ 



To restore the perturbativity in the theory we can take either of the two approaches (i) 
add new particles with new Yukawa contributions so as the keep the Yukawa couplings 
perturbative all the way up to the GUT scale (ii) take the view point that the theory is 
valid only up to a scale which is allowed by perturbative constraints and then some new 
non-perturbative physics takes over. We will study the second option in the present work 
and note that the first option has already been considered by others. To this end, we do 
a complete scan of the allowed regions of the high scale, Mx, and tan/3 for a given set of 
values of m t >,my,m T i. In |Fig. 3 on the next page] we present the results using the lower 
limits on the fourth generation masses given in Eq. (2) and Eq. (3) on page on page HI We 
have demanded that all the Yukawa couplings remain perturbative up to the scale Mx- We 
have presented the results using both 1-loop as well as 2-loop RGE to show the importance 
of using 2-loop RGE. From the 1-loop plot we see that for tan /3 ~ 2, the Yukawa couplings 
barely remain perturbative up to 8 TeV or so. At the two loop level this scale increases to 
about 16 TeV. 

We have also compared our analysis with that of Ref. [20]. We have chosen the same 
masses m t i = 220 GeV, my = 190 GeV and m T > = 100 GeV. The allowed regions in Mx 
vs. tan/3 plane are presented in Fig. 3j|(b) on page on the next page. We find that while 
qualitatively we agree with them, quantitatively we differ in the maximal value of Mx by 
a couple of orders in magnitude. 



on page on the following page shows that if we choose the fourth generation 
fermion masses to be at the weaker lower limits (see Eq. (3) on page 4), we can have 
perturbativity all the way up to ~ 1000 TeV. Note that these weaker lower limits are 
obtained under very reasonable assumptions. Finally, the relevant Mx and tan/3 ranges 
can be read out from the plots as follows : 



1 % tan /3 < 3.6, 1 TeV < M x < 16 TeV for |Eq. (2)| on page S 

1< tan/3 < 2.4, 1 TeV < M x < 920 TeV for|Eq. (3)|on pagel (12) 



3.1 Implications for mSUGRA 

Before closing this section, let us comment on the possibility of realising minimal super- 
gravity with four generations. From a phenomenological point of view such a possibility 
would be interesting with supersymmetric partners of the fourth generations leading to 
new experimental signatures at the weak scale. Further, due to the presence of additional 
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10 12 14 16 18 20 
M x [TeV] 



(a) For the masses in Eq. (1) on page 3 without T-parameter constraints 




max(h M ,h b4 ,h 4 ) at 1-loop and " G = 4 | 
= 220 GeV, m = 190 GeV, m„ = 100.0 GeV ' 



max(h M ,h b4 ,h 4 ) at 2-loop and n Q = 4 

220 GeV, = 190 GeV, m_ 4 = 100.0 GeV 




(b) For the masses considered in Ref. [20]. 




(c) For the masses in Eq. (3) on page 4 with T-parameter constraints 



Figure 3: Regions in the Mx~tan (3 plane where all the Yukawa couplings of the charged fourth generation 
fermions are perturbative at 1-loop and 2-loops respectively. The different rows correspond to different 
values of the fermions indicated on the panels and explained in the subcaption of the figure in each row. 



10 



Yukawa couplings, the weak scale supersymmetric mass spectrum would most likely be 
quite different from that of mSUGRA3 [18]. If the fourth family Yukawa couplings are 
large, they could contribute significantly to the lightest Higgs mass at the 1-loop level, 
thus alleviating the little hierarchy problem [23], [51]. However, in the standard picture 
of mSUGRA with universal or non-universal soft masses at the GUT scale, radiative elec- 
troweak symmetry breaking induced by the large t' (and possibly V , r') Yukawa coupling 
and the neutralino as the lightest supersymmetric particle, would require the theory to 
remain perturbative all the way up to the GUT scale. From the analysis above, we see 
that theory would remain perturbative with four generations only if the fourth generation 
masses are much lower than their experimental lower limits, in fact, closer to the third 
generation masses. In the following we will consider a toy model of mSUGRA4, where the 
masses of the fourth generation are set equal to the their third generation counterparts: 



writ, Tfiy = nib, m T i 



We then compute the mass spectrum at the weak scale for a sample point and compare it 
with that of mSUGRA3. We use our new tool INDISOFT to compute the spectrum at the 
weak scale. 



mSUGRA and n„ = 3 
■> 800 r "V 100 " m = 250 ' A =0 > tan P= 10 > s 9 n 



■ Gauginos 
- Squarks 
. Sleptons 




10 3 10 4 io 5 io 6 io 7 io 8 io 9 ™ 10 ™ 11 ™ 12 ™ 13 ™ 14 ™ 15 ™ 16 ™ 17 ™ 18 

Q [GeV] 

(a) Three generations. 



mSUGRA and n Q = 4 

m^-lOO, m =250, A =0, tan £=10, sgnn=+ 







Gauginos 




Squarks 




Sleptons 




1 — — «Jk hi 2 + r 




10 3 10" 10 s 10* 10 7 1 8 1 g 10 1 °10"l0 12 10 13 10 14 10 1s 10 1 *10"l0 ,a 

Q [GeV] 



(b) Four generations. 



Figure 4: The running of the various soft masses in the MSSM3 and MSSM4 is shown in the left and 
right panel, respectively. The legend is indicated in the panel. The Higgses are the black lines starting at 
V^ 2 + ™o- 



The very presence of four generations, irrespective of whether their Yukawa couplings are 
large or not, can lead to interesting differences between the three and four generation RG 



evolution of the soft masses. In Fig. 4 we show the running patterns of the soft terms 
for three as well as four generations. As expected, the Higgs mass terms (black line) 
run to more negative regions as compared to three generation case, and in fact, they can 
both become negative even for small tan f3 in the four generation case. Perhaps the most 
interesting aspect is the running of the gluino mass, which now due to a smaller /3-function3, 
almost does not evolve (up to 1-loop level) in the four generation case (thick red line). The 



7 Note that h = 3 for MSSM3 and 6 3 = 1 for MSSM4. 
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sleptons (dashed-blue) are not significantly affected, however the squarks (undashed-blue) 
run to lighter values due to reduced gluino running effects. All these differences in the 
running of the soft terms would make themselves evident in the mass spectrum at the 
weak scale. 



Higgses [GeV] 


Gaug 


inos [GeV] 


Squarks & Sleptons [GeV] 


h° 106.7 




96.6 




568.2 


h 


587.4 


A 382.2 


X\ 


178.3 


UR 


547.5 


h 


411.0 


H° 382.6 


xi 


343.0 


d L 


573.6 


bi 


519.9 


H± 390.9 


xl 


362.8 


d R 


546.6 


h 


547.2 




xt 


178.0 


e L 


205.7 


Tl 


209.1 




xt 


364.5 


e R 


146.7 


T 2 


138.9 




g 


607.0 


v e 


189.8 


V T 


189.1 



Table 1: MSSM spectrum with 3 generations and mSUGRA boundary conditions: m = 
100 GeV, mx/2 = 250 GeV, A = GeV, tan/3 = 10, sgn/i=+. The unification scale is 
M GUT = 2.40 x 10 16 GeV. 



Higgses [GeV] 


Gaug 


inos [GeV] 




S quark 


3 & Sleptons [GeV] 


h° 


119.5 


X? 


44.1 


u L 


480.4 




499.7 


t'i 


498.8 


A 


486.5 


X°2 


83.4 


u R 


462.6 


h 


357.8 


t> 2 


356.4 


H° 


486.2 


xl 


474.2 


d~L 


486.7 




432.4 


b[ 


428.7 


H± 


492.8 


xl 


478.1 


d R 


462.0 


h 


465.9 


b> 2 


466.2 






xt 


83.4 


e L 


187.7 


Tl 


196.4 


t[ 


196.2 






xt 


481.4 


e R 


142.0 


T 2 


126.5 




127.1 
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352.1 




170.4 


V T 


169.6 


K 


169.6 



Table 2: MSSM spectrum with 4 generations and mSUGRA boundary conditions: mo = 
100 GeV, mi/2 = 250 GeV, A = GeV, tan/3 = 10, sgn/i=+. For the theory to be 
perturbative, we have chosen all 4th generation masses to be equal to their 3rd generation 
counterparts (toy model). The unification scale is M GUT = 8.82 x 10 16 GeV. 

In ITab. II and ITab. 2\ we present the weak scale spectrum for some sample point. Com- 
paring the spectrum in the two tables, we find the following: (a) the Higgs mass is heavier 
in mSUGRA4, in fact above the LEP limit (b) the lighter neutralinos in mSUGRA4 are 
lighter compared to mSUGRA3, with one close to 60 GeV (c) the squark and the gluino 
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masses are also significantly lighter compared to mSUGRA3 (d) slepton masses do not 
have much of an impact and they seem to be close to those in mSUGRA3. 

Thus the addition of a fourth chiral generation to MSSM3 seems to give rise to a lighter 
super symmetric spectrum at the weak scale with a less fine tuned Higgs mass. These re- 
sults, though valid only in this particular toy model, seem to indicate the possible features 
the supersymmetric spectrum would have if one could make theory perturbative. As dis- 
cussed earlier, one possible way would be to add additional vector-like matter [20]. An 
alternative approach, which we follow, is to lower the scale of supersymmetry breaking and 
ask whether the above features are replicated. We thus look for a low scale supersymme- 
try mediation mechanism which is preferably close to the non-perturbative regime in the 
Yukawa couplings. 

4 GMSB with Four generations 
4.1 Minimal Messenger Model 

Given that the four generational MSSM is barely perturbative up to few hundreds of TeV, 
supersymmetry breaking should be communicated within this energy scale to the visible 
sector. Gauge mediated supersymmetry breaking (for a review, see Refs. [52] and [53]) is 
one such possibility where mediation scales can be as low as 10-20 TeV. In the present 
section, we explore the possibility of reconciling the minimal messenger model (MMM) of 
GMSB with the four generational MSSM. 

The minimal messenger model has a set of chiral superfields (messengers) which transform 
as fundamentals under SU(5). They are coupled to a singlet field S which parametrises 
the hidden sector supersymmetry breaking and whose F-component and scalar component 
attain vevs. As a consequence, messengers attain both supersymmetry conserving and 
supersymmetry breaking masses. This breaking information is then passed on to the 
MSSM sector through gauge interactions. The gauginos attain masses at the 1-loop level, 
given by 

Mi(X) = a l (X)Ag(x) (13) 

where X ~ (S) is the messenger mass scale (up to a coupling constant), A = (F s ) / (S) and 
i — 1,2,3 for the three gauge groups of the MSSM. The tilde on the gauge couplings a; 
denotes a division by 4ir. g(x) is the loop function with x = (S) 2 / (Fs). This function rises 
monotonically between g(0) = 1 and g(l) = 1.386. We chose x = 1/2 where g(x) ~ 1. 
The scalars attain their masses at the two loop level and they are given by : 

3 

mj(X) =A 2 J^Ci d 2 (X) h(x) (14) 
i=i 

where X is the messenger scale, C 3 = 4/3 for SU(3) triplets and C 2 = 3/4 for SU(2) 
doublets and C\ = Y 2 represents the hypercharge of the particles. The function h(x) 
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which is almost flat for most x values is also equal to 1 at x — 1/2. / runs over all the 
scalars in the theory, / = {Q, U, D, L, E, Hi, H 2 }. Going from three generations to four 
would not change any of these boundary conditions at the messenger scale X. However, the 
weak scale spectrum is expected to be completely different from that of the three generation 
case due to the presence of large fourth generation Yukawa couplings. It is interesting to 
study the soft spectrum in the presence of large Yukawa couplings, but a small running 
scale ~ log (Mx/M SVSY ) « 4 — 5. To calculate the low energy spectrum we need to compute 
the RG effects on the soft terms and evaluate the soft mass matrices at the weak scale. To 
that end, we use INDISOFT which is described in more detail in Appendix [B] 
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Figure 5: Regions in MMM parameter space A-M mess . The lower-diagonal part is ruled out as A > M mess . 
In the upper-diagonal part, from left to right, the first region (red) tachyonic r', and the second (orange), 
third (cyan), fourth (green) do not have consistent radiative electroweak symmetry breaking as indicated 
by the tachyonic Higgses. 



Choosing the masses for the fourth generation to be given by Eq. (3) on page 4 and tan j3 = 
1.75, we have computed the supersymmetric spectrum at the weak scale. We vary 10 TeV < 
A, M mess < 500 TeV. The results are characterised by various regions in the parameter 
space as shown in |Fig. 5[ The blue region in the lower diagonal part is ruled out, because 
the messenger scale M mess is smaller than A. The greater part of the parameter space does 
not have radiative electroweak symmetry breaking as indicated by the tachyonic Higgs 
scalars (orange, cyan, and green regions) and is thus ruled out. There is, however, a small 
red region where electroweak symmetry breaking is possible, but the r' is a tachyon. 

It is interesting to study in detail the region of the tachyonic r' , with 
10 TeV < A < 70 TeV, 30 TeV < M x < 300 TeV. 



A sample spectrum of this region is given in Tab. 3 on the next page where we have chosen 
A = 50 TeV and M mess = 100 TeV. Notice that the light Higgs mass does not satisfy the 
LEP-II constraints for this point in the parameter space. 
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Higgses [GeV] 


Gaug 


inos [GeV] 




S quark 


3 & Sleptons [GeV] 


h° 


46.2 


v9 

Al 


64.3 


Ut 


758.1 




766.1 


t\ 
h \ 


722.6 


A 


507.6 


A2 


127.0 


^ it 


735.5 


h 


639.3 


t> 2 


583.8 


H° 


532.2 




640.6 


d L 


761.1 


\ 


725.1 


b[ 


733.4 


H± 


516.1 


xl 


655.1 


d R 


733.8 


h 


734.3 


b> 2 


525.5 








126.9 


e L 


208.3 


T\ 


208.4 


r[ 


320.3 








652.0 


e R 


88.1 




87.8 


% 


193.4 






<? 


438.4 




197.2 


v T 


197.2 




202.7 



Table 3: Minimal GMSB spectrum with 4 generations: — 1 GeV, M mess = 100 TeV, 
A = 50 TeV, tan/5 = 1.75, sgn/i=+. f is tachyonic, m h = 46.2 GeV, mg = 1.2 x 10~ 9 
GeV (gravitino), NLSP is neutralino. 



To trace the reasons for a tachyonic r', let us consider its mass matrix which is given by 



m~, 



m 
m 



LL 

2 * 
LR 



III 

m 



LR 
2 

RR 



o 



ml 
ml |U ' 

T 2 



(15) 



where 



/ 1 

m LL = m~,^ + m T t — m z cos2p I - — sin &w J ~ m?^ 



= + m T , — m z cos 2/3 sin fc'iy ~ m 7 i 

m 2 LR = v d (A*, — jjt, y T > tan j3) ~ — m r //x tan (3 (16) 



where O is some mixing matrix which diagonalises the r' mass matrix. Unless both r' are 
tachyonic, the condition that there is no tachyon in the spectrum is given by: 



m~^m~i — (m T i /itan/3) 2 > (17) 



In the MMM model, it is typical that fi is very large, between one TeV to tens of TeV. In our 
analysis, [i is determined at the weak scale by electroweak symmetry breaking conditions. 



With m T t = 75 GeV, the second term of Eq. (17) is quite large. After the RG evolution, 
the soft terms for Tl,r take the form [54]: 



m 2 ~,(qZ) « m 2 l (X)-a T/ (2mf(X)+m^(X))\og(X/qZ) 



mUqZ) « m 2 ~(X)-2a T/ (2m 2 T (X)+m 2 ~(X))\og(X/qZ) 



where ay = y 2 ,/ (167r 2 ) « 1.1 x 10 3 x (1 +tan 2 /9), qZ « 1 TeV is the typical supersym- 
metric mass scale close to the weak scale, and X represents the messenger scale as before. 
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Note that the log-factor is extremely small, ~ 4.6. Thus, the effective difference between 
the high scale and the weak scale sleptons is very small, ~ (5 — 10)%. For this reason, 
it becomes difficult to maintain Eq. (17) on the previous page positive. Thus, the MMM 
model is perhaps too restrictive for the case of four generations. With modified boundary 
conditions at the messenger scale, it might be possible that gauge mediated supersymmetry 
breaking could lead to a phenomenologically viable spectrum. This brings us to the new 
ideas of general gauge mediation (GGM). 



4.2 General Gauge Mediation 

In the last couple of years, there has been a paradigm shift in the way we understand gauge 
mediated supersymmetry breaking. Ref. [55] has introduced the framework of general gauge 
mediation where the starting point is that the hidden sector and the visible MSSM sector 
should completely decouple in the limit where the gauge couplings are set to zero. The 
generalisation lies in the set of formulae for the soft terms which span from models with 
a weakly coupled messenger sector as well as to the strongly coupled ones. The general 
gauge mediation can allow for different scales for the scalars and the gaugino masses, with 
As f° r scalar masses and for gaugino masses. In fact, the authors of Ref. [56] argue 
that with this kind of parametrisation, the GMSB analysis is closer to mSUGRA where 
the role of Aq is played by m 1( / 2 and that of As is played by m . Such freedom in the 
choice of scales could lead to a relaxation of the parameter space. Finally, let us note that 
this framework of general gauge mediation does not have any solution to the /i-problem 
as it is not generated through the above interactions. Typical solutions of the /i-problem 
could decouple the Higgs sector from the sleptonic sector which might pave way for the 
solution to the tachyonic r problem. A more detailed exploration regarding these issues 
will be addressed in a forthcoming publication. 



5 Summary and Outlook 

In the present work, we have revisited in detail the perturbativity constraints on the 
fourth generation Yukawa couplings in the minimal supersymmetric standard model. We 
have shown that if the fourth generation masses lie close to their present upper limits 
(or within 25% of these limits), it is possible to have perturbativity only up to ~ 1000 
TeV; tan/3 is confined to be very low in these models, ~ 3. This makes it hard to reconcile 
gravity mediated supersymmetry breaking models with four generations unless one assumes 
additional vector-like matter. However, as demonstrated in the spectrum of the toy model 
of mSUGRA with four generations, the presence of a fourth generation can lead to new 
features like lighter gluino and squark masses which could be accessible at LHC. 

With this motivation, we have studied gauge mediated supersymmetry breaking with four 
generations. We have explored the already highly constrained minimal messenger model 
with the ranges allowed and found that in most of the parameter space, radiative elec- 
troweak symmetry breaking is not possible. There is a small region where electroweak 
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symmetry breaking is viable but r' is tachyonic. It is difficult to accommodate a fourth 
generation in the highly constrained minimal messenger model. We speculate that the more 
general framework of general gauge mediation would be suited to get phenomenologically 
viable mass spectra with four generations. 

The combination of supersymmetry and four generations seems to lead to interesting su- 
persymmetric spectroscopy at the weak scale. In the present work, we have just explored 
possible supersymmetry breaking scenarios which can work within this framework. How- 
ever, four generational MSSM would have strong implications in various other sectors like 
flavour physics, neutrino physics and most interestingly the Higgs physics in MSSM. These 
aspects as well as concrete models within gauge mediation need to be studied. Experi- 
mental consequences of the presence of the fourth generation with SUSY, in B-physics, 
D-physics as well as direct detection at LHC needs to be investigated. 
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A The Renormalisation Group Equations 



For completeness, we list the 2-loop renormalisation group equations for the gauge and 
Yukawa couplings that we have used to generalise SOFTSUSY to the case of the MSSM with 
hq chiral generations. These RGEs can be derived from the general case [57-59] and are 
well-known in the literature [16,18,60-64]. Here we are following Ref. [64]. 



A.l The Gauge Couplings 

The running of the gauge couplings is given by 



dt 



16tt 2 



bd + 



(16tt 



2^2 



hk9l - Tr {c iu Yl Y u + C id Y d Y d + C ie Y\Y e ) 



, (19) 



where t = log(/x//xo), i s the renormalisation scale, fj,Q some reference scale, and i,k = 
1,2,3 refer to the gauge groups U(l)y, SU(2)l and SU(3) C . At 1-loop, the 6, are given by 



~ " 2n G , 
5 



2n ( 



9 - 2n G . 



(20) 



The fourth generation comes in complete GUT multiplets, and thus most predictions from 
grand unification are unchanged. In fact, the only modification at 1-loop is that a GVT 
increases. At 2-loop, the Yukawa couplings enter the RGEs, and the coefficients are given 
by 



(b) 



ik 



/ 38 


6 


88 \ 






5 


15 1 




I 


14 


8 


H 




3 


68 J 




\ 15 




3 / 





9 

f 
5 


9 
5 

-17 












-54 



26 14 18 

5 5 5 

6 6 2 
4 4 



(21) 



where i, k — 1, 2, 3 as before, and f = u,d, e. 



A. 2 The Yukawa Couplings 

The running of the Yukawa couplings is given by 

/ 1 „m 1 



-Y = Y I - ft^ 4- ft^ 1 (22) 

^ *- u,d,e x u,d,e I 16lT 2 ' d ' e (167T 2 ) 2 u ' d ' e ' 



At 1-loop, the /^-function coefficients are: 



flf' = 3 V i V. + YlV < + 31> { YiY.}-g 9! + M + ^) 
/#> = 3YjY d + YtY„ + Tr {syJy, + YtY e } - (lg* + 3^ + fg 2 3 
(3^ = 3YtY e + Tr {syJy, + YjY e } - + 3^ 2 



(23) 
(24) 
(25) 
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Note that the number of generations enters only through the dimensionality of the Yukawa 
matrices. At 2- loop, tlq is explicitly present, and the /3-function coefficients are: 

/9?> = -4(YtYj 2 -2(YjY d ) 2 -2YjY d YtY u -9Tr{YtYjYtY u 

- Tr {sY^Y, + YtY e } Y^Y, - 3Tr {3(Yt Y J 2 + YjY d YtY u } 

+ {\& + QglM^u + (^i 2 )YjY, + {~gl + 16^ 2 )Tr {Y+Yj 



,26 403 4 21 4 32 304, 

T5 nG + lEo )91 + {6nG ~ + { Y nG ~ — 



+ ghl + ^ght + 8 9 W 3 (26) 

/3f = -4(YtY d ) 2 -2(YtYj 2 -2YtY„YtY (i -3Tr{YtY. il }YtY u 

- 3Tr {3YjY d + YlY e } YjY d - 3Tr {s(yJy,) 2 + (Y^YJ 2 + Y^Y^Y, 



+ (^i)YlY u + (~g 2 + e^YjY,, + (--g 2 + 16^ 2 )Tr {yJY,} + (-^ 2 )Tr {Y|Y e } 
+ {^n G + ^)<? 4 + (6n G - y )^ 2 4 + (yn G - + g\g 2 + + 8g 2 2 g 2 3 (27) 

0?> = -4(YtY e ) 2 -3Tr{3Y;Y d + YtY e }YtY e -3Tr{3(YjY,) 2 + (YtY e ) 2 

+YjY d Yt Y u ] + (6^ 2 )YtY e + (^ 2 )Tr {Y e f Y e } + (-^ 2 + 16^ 2 )Tr {yJy,} 

+ (fn G + 2 l)gi + (6n G - |)^ + jj^ (28) 

A. 3 Approximate Upper Bounds 

An approximate upper bound on the fourth generation t' can be obtained in the limit where 
m t i dominates over all the other masses, m t i 3> {m t , my, m T /} purely from perturbative 
requirements. To see this, let us rewrite Eq. (22) on the preceding page in terms of the 
dominant Yukawa couplings: 

K = h'Uhl + hl + 3h 2 -^gi-3gl-^g 2 

Phy = h v Lh 2 +h 2 +hl,-^gl-3gl-l-g 2 



K = h T ,(3h 2 +4h 2 T ,-3g 2 2 -^g 2 ^ 



Pht = h\3h 2 , + Zh 2 -^g 2 -3g 2 -^g 2 \ (29) 

These equations are non-linear and coupled and obviously, should be solved numerically. 
However an analytical estimate can be obtained in the above mentioned limit m t i ^> 
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{m t , my, m T r}. In this limit the solution for hp can be written as [65] 



kf(z) 



E(z)k tl (0) 



:i + Qk tl (0)F(z)Y 



(30) 



where z = 2log(Mx/Mz) with Mx denoting the high scale and Mz, the weak scale. 
kp = h\, j (47r) and 



E(z) 
F(z) 



13 S_ JJ5_ 

, ~ ttl(0) \ 9»1 / r Qt 2 (0) \ 6 2 / ~ CK3(0) \ 36 3 

I + l + 6 2 — ^z) 1 + 63-^2 1 



47T 



47T 



47T 



(31) 



fcf(0) is the value of the t' Yukawa (squared) at the high scale, Mx- These are the same 
expressions one obtains for the top quark Yukawa coupling within three generations. In 



four generations, the only difference is the /3-functions of the gauge couplings. In Eq. (31) 

43/3. If the Yukawa 



— £>i, and from Eq. (20) on page 18, b 



-1, bo 



3, 61 



coupling becomes very large, from Eq. (30) , we can derive an upper bound on the (square) 
of the Yukawa of the top-prime in the limit k t >(0) — > 00: 

E(z) 



k t '(z) 



6 F(z) 



(32) 



In lTab. 41 we present the upper limits on the m t i in this approximate limit using Eq. (32) 
We see that the present limit on the m t > readily rules out perturbativity up to the Planck 
scale or even the GUT scale. If we demand perturbativity beyond 1-100 TeV, mp is forced 
to be close to its experimental present lower bound or below. Beyond 10 6 GeV the present 
limit already rules out perturbative a Yukawa coupling for the t'. As we have seen in the 
main text, demanding perturbativity of the b' and r' Yukawa couplings would put further 
stringent constraints on the scale Mx- 



M x (GeV) 


10 18 


10 16 


10 11 


10 6 


10 3 


m t , (GeV) 


211.3 


212.8 


222.2 


264.2 


437.7 



Table 4: The approximate upper limits on mp in the limit where mp dominates over all 
the masses, for various high scales, Mx- 



B INDIS0FT 

In this appendix, we summarise the most important features of the software tool that we 
have developed to calculate soft spectra in the MSSM with four generations. More details 
will be presented in a separate publication. 
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INDISOFT is based on SOFTSUSY 3.0.9. [50]. We have preserved the original C++ class struc- 
ture and extended the functionality of the individual classes that handle the calculations. 
It is beyond the scope of the present appendix to describe SOFTSUSY in detail, and we refer 
for that to its manual [50]. 

The fermion masses and gauge coupling constants are stored in the class QedQcd which we 
have extended to comprise also the fourth generation fermion masses. The light fermion 
masses and the couplings that are entered at different scales (depending on where they are 
experimentally known) are then run to Mz using the class RGE that implements the renor- 
malisation group running for every class that derives from it, in this case for QedQcd. For 
the heavy fermions, we have generalised the functions that calculate the pole mass from the 
running one and vice versa. From now on we will assume that we are working in a frame- 
work (like mSUGRA or gauge mediation) where the soft masses are obtained from boundary 
conditions set at a higher scale. The class MssmSusy is derived from RGE and contains the 
supersymmetric parameters of the theory; the Yukawa matrices have been promoted to 
4x4 matrices to include the fourth generation. The /^-functions of MssmSusy have been 
generalised to the case of four generations (see Appendix [Aj. The class SoftParsMssm is 
derived from MssmSusy and contains the soft supersymmetry-breaking terms of the MSSM; 
the soft mass matrices and the trilinear couplings have been generalised; the supersym- 
metric parameters and their RGE evolution are inherited from MssmSusy; the /^-functions 
of SoftParsMssm for the soft terms have been generalised. The class MssmSof tsusy derives 
from SoftParsMssm and organises the actual calculation of the spectrum. QedQcd is used to 
initialise/guess the SUSY parameters at the scale m t which are then run by RGE to the high 
scale Mx where the boundary conditions on the soft terms are imposed. Note that for the 
gauge and Yukawa couplings, we use 2- loop RGEs (see Appendix |A|, and 1-loop RGEs 
for the rest. MssmSof tsusy is then run back to M SUSY (= Mz for the first iteration), where 
the electroweak symmetry breaking conditions are checked and the physical spectrum is 
calculated (at tree-level for the first iteration, and later at 1-loop). The class sPhysical 
that calculates the physical masses and stores the results has been generalised to handle 
four generations. The class drBarPars that derives from sPhysical and manages the DR 
parameters for the calculation has been generalised. The tadpoles, 1-loop radiative correc- 
tions to the Higgses, squarks, sleptons, and the threshold corrections to the gauge couplings 
have been generalised to include contributions from the fourth generation fermions. The 
aforementioned steps are iterated until satisfactory convergence is achieved. 

As a result of our work, we found some minor typos and bugg^l in SOFTSUSY that have been 
fixed in subsequent versions. We have rewritten the linear algebra classes from scratch^], 
replacing the pointer constructions used to represent vectors/matrices and the algorithms 
operating on them by the Standard Template Library (STL) containers and algorithms. In 

8 S0FTSUSY 3.0 to 3.0.9 did not correctly indicate the regions where no electroweak symmetry breaking 
is possible, and in release 3.0.9, there were some typos in the formulae for the radiative corrections that 
affected the calculation of the physical masses at the subpercentage level. 

9 During the final stages of this publication, we became aware of the SOFTSUSY 3.1 release in which the 
linear algebra classes have been rewritten by D. Grellscheid. We have not compared our changes to his. 
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addition to the presently available mSUGRA and minimal GMSB boundary conditions, 
we have implemented a model of general gauge mediation along the lines of Ref . [56] . We 
have linked ROOT [66] to our programs to generate plots both interactively and in batch 
mode. In future, we plan to extend INDISOFT by right-handed neutrinos. Due to space 
limitations, we must refrain from discussing all the changes we have made. 
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